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Introduction 


The  Cauchy  stress  tensor  T  can  be  decomposed  into  a  spherical  part,1 
—pi,  and  a  deviatoric  part,  T*: 

T  =  — pl  +  T*  ,  ti  =  -p  +  t?.  (1.1) 

The  pressure  p  is  given  by 

p=-ItrT  =  -|(<1+t2  +  i3),  (1.2) 

and  the  deviatoric  stress  tensor  T*  is  given  by 

T*  =  T  -  |(tr  T)I ,  tt  =  U  -  !(<!  +t2  + 13) ,  (1.3) 

so  that 

tr  T*  =  t{  + 12  +t3  =  0 .  (1.4) 

The  principal  stresses  t-  and  principal  deviatoric  stresses  t*  axe  the  princi¬ 
pal  values  of  T  and  T* ,  respectively.  For  an  isotropic  elastic  material  the 
linear  theory  predicts  that  the  pressure  depends  only  on  the  volumetric 
strain,  whereas  T*,  which  is  a  tensorial  measure  of  shear  stress,  depends 
only  on  the  shear  strain. 

The  object  of  this  report  is  to  study  those  aspects  of  the  nonlinear 
elastic  response  of  polycrystalline  metals  and  ceramics  which  may  be  im¬ 
portant  in  high  velocity  impacts,  where  large  elastic  increases  in  density 
are  encountered.  As  is  well-known,  a  nonlinear  pressure/density  relation 
is  generally  required  for  such  applications.  Here  we  focus  on  other  nonlin¬ 
ear  effects  that  are  often  ignored  or  treated  incorrectly.  In  particular,  we 
study  the  coupling  of  the  pressure  and  deviatoric  stress  which  results  from 
the  dependence  of  shear  stress  on  volumetric  strain  and  the  dependence  of 
pressure  on  shear  strain. 

We  begin  with  a  brief  review  of  relevant  results  from  the  theory  of  hy¬ 
perelastic  materials.2  Let  F  denote  the  deformation  gradient  relative  to 
some  fixed  reference  configuration.  An  elastic  (or  Cauchy-elastic)  material 
is  one  for  which  T  is  a  function  of  F  only.  A  hyperelastic  (or  Green-elastic ) 
material  is  an  elastic  material  for  which  the  first  Piola-Kirchhoff  stress  ten¬ 
sor  (det  F)T(F-1)T  is  the  gradient  of  some  real- valued  function  e  of  F,  in 

1  Second-order  tensors  are  denoted  by  boldface  uppercase  Roman  letters.  The  identity 
tensor  is  denoted  by  I,  and  tr  denotes  the  trace  function.  The  deviatoric  part  and  the 
transpose  of  any  tensor  A  are  denoted  by  A*  and  A1",  respectively.  The  norm  of  A  is 
II  A||  =  tr  (AtA). 

2 A  standard  reference  is  Truesdell  and  Noll  [1], 
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which  case  e  is  called  the  strain  energy  or  the  stored  energy  per  unit  refer¬ 
ence  volume.3  We  consider  only  hyperelastic  materials  here.  For  simplicity 
we  also  assume  the  material  is  isotropic.  These  conditions,  together  with 
the  requirement  that  constitutive  equations  be  properly  invariant  under 
changes  of  the  frame  of  reference,  yield  various  reduced  forms  (see  below) 
for  the  strain  energy  and  the  Cauchy  stress  tensor. 

The  left  stretch  tensor  V  is  the  unique  symmetric  positive-definite  ten¬ 
sor  occurring  in  the  left  polar  decomposition  of  the  deformation  gradient: 

F  =  VR,  V2  =  B  =  FFr ,  (1.5) 

where  R,  the  rotation  tensor ,  is  proper  orthogonal,  and  B,  the  left  Cauchy- 
Green  tensor ,  is  symmetric  positive-definite.  The  tensors  V  and  B  share 
a  common  set  of  principal  axes,  called  the  principal  axes  of  strain  in  the 
deformed  state  or  the  Eulerian  strain  axes.  The  principal  stretches  A,- 
(i  =  1,2,3)  are  the  principal  values  of  V,  and  the  principal  strains  are 
A,-  —  1.  The  principal  values  of  B  are  &,•  =  A?.  Let  J  denote  the  Jacobian 
of  the  deformation,  and  let  p  denote  the  ratio  of  the  densities  p  and  p0  in 
the  deformed  and  reference  configurations,  respectively: 


J  =  detF  =  det  V  =  AXA2A3 ,  p  =  p/p0  =  l/J.  (1.6) 


Relative  to  an  undistorted  reference  configuration  of  an  isotropic  hypere¬ 
lastic  material,  e  and  T  may  be  expressed  as  isotropic  functions  of  V  or 
B;  and  the  Eulerian  strain  axes  are  principal  axes  for  T.  e  may  also  be 
expressed  as  a  symmetric  function  of  Ax,  A2,  A3  or  bu  b2,  b3.  Then4 


de  de 

T^V-  =  2^B_, 


de  de 

(1.7) 


Note  that  the  summation  convention  is  not  used  here  or  elsewhere. 

3This  restriction  on  the  constitutive  equation  is  motivated  by  the  restrictions  imposed 
on  thermoelastic  materials  by  the  second  law  of  thermodynamics.  A  thermoelastic  ma¬ 
terial  reduces  to  a  hyperelastic  material  for  isentropic  deformations,  provided  we  take 
e  to  be  the  internal  energy  per  unit  reference  volume.  The  same  thermoelastic  mate¬ 
rial  reduces  to  a  (generally  different)  hyperelastic  material  for  isothermal  deformations, 
provided  take  e  to  be  the  free  energy  per  unit  reference  volume. 

4If  e  is  a  scalar- valued  function  of  a  symmetric  tensor  A,  then  the  gradient  of  e  at  A 
is  the  symmetric  tensor  de/d  A  with  the  property  tr  [(de/<9A)S]  =  d/dt  e(A  +  *S)|t=o, 
where  S  is  any  tensor  such  that  A  +  tS  lies  in  the  domain  of  e  for  sufficiently  small  t. 
Furthermore,  e  is  an  isotropic  function  of  the  symmetric  tensor  A  iff  e  is  a  symmetric 
function  of  the  principal  values  a,  of  A,  in  which  case  de/dA  is  coaxial  with  A  and  has 
corresponding  principal  values  de/dai . 
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2  The  Logarithmic  Strain  Tensor 

The  principal  logarithmic  strains  li  are  the  logarithms  of  the  principal 
stretches,  and  the  (Eulerian)  logarithmic  strain  tensor  L  is  the  tensor 
coaxial  with  V  and  B  but  with  corresponding  principal  values  equal  to 
the  principal  logarithmic  strains: 


/,•  =  In  A;  =  |  In  bi ,  L  =  lnV  =  flnB.  (2.1) 


Then  the  strain  energy  function  e  of  an  an  isotropic  hyperelastic  material 
can  be  expressed  as  an  isotropic  function  of  L  or,  equivalently,  as  a  sym¬ 
metric  function  of  the  principal  logarithmic  strains.  From  (1.7)3  and  (2.1)! 
it  follows  that 


de 

t‘=pai ;■ 


_  „de 

t  =  psl- 


(2.2) 


The  logarithmic  strain  tensor  has  several  interesting  properties.  By 
(2.1)  and  (1.6)  we  have 


tr  L  —  "I-  ^2  "b  ^3  —  ln( Aj  A2  A3)  —  In  3  —  —  In  p .  (2.3) 

Hence  trL  is  a  measure  of  the  volumetric  strain.  The  deviatoric  part  L* 
of  L  has  principal  values 


where  i,j,  k  always  denotes  a  permutation  of  1,2, 3.  Now  consider  another 
deformation  which  differs  from  the  given  one  by  a  superimposed  dilatation 
and  rotation,  so  that  this  second  deformation  has  left  stretch  tensor  aV 
and  principal  stretches  a\1,a\2,a\3  for  some  a  >  0.  Then  L*  is  the  same 
for  both  deformations,  since  the  factor  a  cancels  out  in  (2.4).  Hence  the 
deviatoric  logarithmic  strain  tensor  L*  is  independent  of  the  dilatational 
part  of  the  deformation  gradient,  i.e.,  it  is  independent  of  the  volumetric 
stretch  or  strain  as  measured  by  (1.6)  or  (2.3).  Therefore  L*  is  a  tensorial 
measure  of  shear  strain  only. 

From  the  results  above  it  follows  that  e,  T,  T* ,  and  p  may  be  regarded 
as  functions  of  the  independent  variables  L*  and  J  (or  p),  with  each  func¬ 
tion  isotropic  in  L*  for  fixed  J  (or  p).  Then  it  can  be  shown  that  (2.2)2  is 
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equivalent  to  the  conditions5 

» =-(§)L.  -d  jt‘=(^)/  <2-5> 

where  the  subscript  denotes  the  variable  held  constant  during  differentia¬ 
tion.  Equivalently, 

'-'■(SI  *■-'©,•  (2-6) 

In  terms  of  the  pressure  response  functions 

p  =  p(J,V)=P(p,V),  (2.7) 

equations  (2.5)j  and  (2.6)!  may  be  integrated  to  yield 

e  =  £  p((,  L  ’)  df  +  i(L’)  =  L*)e-»  d(  +  i(L')  (2.8) 

for  some  isotropic  function  i  of  L*  only.  From  (2.8)  we  see  that  the  pressure 
response  function  of  an  isotropic  hyperelastic  material  determines  the  strain 
energy  function  to  within  an  arbitrary  function  of  L*  only.  Also  note  that 
for  an  isotropic  elastic  material,  the  existence  of  a  strain  energy  function 
places  no  restrictions  on  the  pressure  response  function.  That  is,  given  any 
isotropic  function  p  of  J  and  L* ,  there  exist  infinitely  many  strain  energy 
functions  e,  namely  those  given  by  (2.8)1 ,  for  which  the  corresponding 
pressure  response  function  is  p  =  p(J,  L*);  then  the  response  function  for 
the  deviatoric  stress  tensor  is  determined  from  (2.5)2. 

In  practice  we  are  usually  faced  with  a  different  situation — namely,  we 
have  a  limited  amount  of  experimental  data  on  the  material  response,  and 
these  data  impose  restrictions  on  the  response  functions  for  both  the  pres¬ 
sure  and  the  deviatoric  stress.  We  then  need  to  determine  any  additional 
restrictions  imposed  on  these  response  functions  due  to  the  existence  of  a 
strain  energy  function.  We  address  this  problem  in  the  next  section  for 
small  shear  strains.  The  following  result,  which  requires  no  restrictions  on 
the  magnitude  of  the  shear  strain,  is  an  immediate  consequence  of  (2.5)  or 
(2-6). _ 

5To  prove  (2.5)i,  for  example,  note  that  since  /*  is  independent  of  J,  (2.4)  yields 
dU/dJ  =  1/(3 J).  Then 


But  from  (1.2),  (2.2),  and  p  =  1/J, 
equal  to  —  p. 


de  dk 
dhdJ 


1  ds 
3 Ifridk' 


we  see  that  the  expression  on  the  right  above  is 
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For  an  isotropic  hyperelastic  material,  the  following  three  conditions  are 
equivalent: 

(i)  The  strain  energy  decouples  additively  into  a  function  of  density 
only  and  a  function  of  shear  strain  only,  i.e.,  £  =  i(p)  -f  e(L*); 

(ii)  The  pressure  depends  only  on  the  density,  i.e.,  p  =  p(p)  ; 

(iii)  JT*  is  independent  of  the  density. 

If  these  conditions  hold,  then  p{p)  =  p2de/dp  and  T*  =  pdi/dL*.  In 
particular,  either  (i)  or  (ii)  implies  T*  is  proportional  to  the  density  ratio 
p.  In  Section  §4  we  show  that  these  properties  axe  generally  inconsistent 
with  experimental  data. 

In  the  remainder  of  this  section  we  consider  some  relations  involving 
the  invariants  of  L*.  These  results  are  utilized  in  the  next  section.  Let  7* 
and  8*  denote  the  second  and  third  moments  of  L*: 

7*  3  tr[(L7M|LT  =  (/,7  +  W)2  +  («! 

- 2 <;,%•  +  i;i;  +  ;,•/,*)  =  2[(/?)»  +  i-i-  +  (//)*] 

=  lift*  -  ft*)2  +  (ft*  -  ft*)2  +  ft*  -  ft*)2] , 

and 

S'  =  tr  [(L*)3]  =  (ft*)*  +  ft")3  +  ft*)3 
=  3detL*=3  l,*W  =  -E(W, 

m^n 

where  we  have  used  the  fact  that 

tr  L*  =  If  +  If  +  I3  =  0 . 

Both  7*  and  8 *  are  isotropic  scalar  measures  of  shear  strain  which  are  inde¬ 
pendent  of  the  dilatational  part  of  the  deformation.  Alternate  expressions 
for  7*  follow  from  (2.9)6  and  the  identities 

I;  -  If  =  ft  -  ft  =  In  .  (2.12) 

Note  that  7*  >  0  and  that  yfy*  =  ||L*||  is  the  norm  of  L*,  which  may 
be  interpreted  as  the  equivalent  shear  strain ;  it  is  zero  iff  the  deformation 
gradient  is  a  dilatation  superimposed  on  a  rotation: 

7*  =  0  &  ||L*||  =  0<S>L’  =  O^V  =  J1/3!  &  F  =  R.  (2.13) 


(2.9) 


(2.10) 


(2.11) 
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Of  course,  7*  =  0  6*  =  0,  but  the  converse  does  not  hold.6 

Let  L**  denote  the  deviatoric  part  of  the  square  of  the  deviatoric  part 
of  L: 

L**  =  [(L*)2]*  =  (L*)2  —  |7*I.  (2.14) 

Then  the  principal  values  1**  of  L**  are  given  by 

ir  =  m2  -  h'  =  jiff)2  +  2(/4*]  =  »;v  +  jr  •  (2.15) 


From  (2.3),  (2.4)1}  (2.9)— (2.11),  and  (2.15),  we  find  that7 
dp  _  dj*  _  55* 


5/,- 


d~f* 

—  =  2/* 

a/i  1 5 


dh 


=  31* 


(2.16) 


As  observed  previously,  any  function  of  L  can  be  regarded  as  a  function 
of  p  and  L*.  Therefore  any  scalar- valued  isotropic  function  of  L  can  be 
expressed  as  a  function  of  p  and  the  second  and  third  moments  of  L*.  In 
particular,  this  applies  to  the  strain  energy  e  and  the  pressure  p: 

£  =  e(p,7*,5*),  P  =  p(P,Y,f)  =  P2  (2.17) 

where  (2.6)j  has  been  used.  Then  from  (1.1),  (2.2)ls  (2.17)3,  and  (2.16), 
we  obtain  the  following  formula  for  the  principal  deviatoric  stresses  in  an 
isotropic  hyperelastic  material: 


equivalently, 


(2.18) 


(2.19) 


Note  that  the  coefficients  of  L*  and  L**  in  (2.19)  generally  depend  on  all 
three  of  the  variables  p,  7*,  and  S*. 

6If  the  principal  stretches  are  ordered  so  that  Ai  <  A2  <  A3,  then 

s*  =  0  h*  =  0  h*  =  -i3m  o/2  =  |(/i  +  h)  o  a2  =  v/aTa7. 

In  particular,  6*  =  0  for  a  simple  shear  (where  Aj  =  I/A3  and  A2  =  1)  and  also  for  a 
simple  shear  with  superimposed  dilatation  (where  Ai  =  J2!3/ A3  and  A2  =  J1/3). 

7In  deriving  (2.16)2,3  it  is  useful  to  first  establish  the  following  results:  dl*  /dh  =  §, 
dl'/dlj  =  -A  dir/dh  =  §/,*,  and  dlT/dlj  =  §/*'. 
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3  Approximate  Formulas 

In  this  section  we  derive  approximate  formulas  for  the  strain  energy,  pres¬ 
sure,  and  deviatoric  stress  that  are  valid  for  arbitrarily  large  changes  in 
density  and  sufficiently  small  shear  strains.  From  (2.17)!,  (2.9),  and  (2.10), 
it  follows  that 


«  =  i(p)  +  e2(p)r  +  Ot(\\V\\3) 

=  e(p)  +  £2(^)7*  +  ez(p) &  +Op{\\L  ||4)  (3.1) 

for  some  functions  e,  e2,  and  e3  of  p  only.  Here  (||L*  ||n)  denotes  a 
function  of  p  and  L*  which  is  of  order  n  in  L*;  i.e.,  there  is  a  function 
M(p)  such  that  C^(||L*||n)  <  M(p)  ||L*||n  for  L*  sufficiently  close  to  0. 
Similarly,  from  (2.17)2  we  have 


P  =  p(p)  +  0A(  ||L*|H 

=  P(p)  +Pj(/5)7-  +  Of  (||L*  ||3) 

=  pW+ftWr+PsWS’+OjdlL-ll^)  (3.2) 


for  some  functions  p,  p2,  and  p3  of  p  only.  Since  7*  and  6*  are  independent 
of  p,  (2.17)3  an<i  (3-1)  imply 


P(P)  =  P2J~'  Pzip)  =  P2-Jjr ,  P3(p)  =  P2^r-  (3.3) 


We  call  p  =  p(p)  the  hydrostatic  pressure  at  the  density  ratio  p.  It  is 
the  pressure  the  material  would  experience  if  the  shear  strain  were  zero; 
i.e.,  if  the  conditions  in  (2.13)  hold,  in  which  case  (1.1)  and  (2.19)  require 
that  T  =  — pi  =  —  p(p)I.  In  accordance  with  experimental  data  for  most 
materials,  we  assume  dp/ dp  >  0.  Then  the  bulk  modulus  k  is  positive: 


0< 


K  = 


(3.4) 


and  p  may  be  regarded  as  a  function  of  the  hydrostatic  pressure  p.  Hence 
any  function  of  />,  such  as  k  and  the  functions  p  and  u>  below,  may  also  be 
regarded  as  a  function  of  p. 

From  (2.18)  and  (3.1),  the  principal  deviatoric  stresses  satisfy 


<•  =  2/.(jS)fi*  +  «%(l|L*|P) 

=  2/,(p)ir+3a,(p)/r  +  Op(||L'||3),  (3.5) 
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where  fx(p),  the  shear  modulus  at  density  ratio  /?,  is  given  by 


and 


M)  =  Mp)  =  P 

^(p)=ps3(p)  =  P  (j^r) 

\  /p,7* 


7*  =6 =0 


7*=5*=0 


The  relations  (3.5)  are  equivalent  to 


T'  =  2^p)V  +  Op\\Vf) 

=  2^)r  +  3^)L"  +  <a(||L-||3). 


(3.6) 

(3.7) 


(3.8) 


From  the  constitutive  relations  (3.2)  and  (3.5)  or  (3.8),  we  see  that  in 
addition  to  the  dependence  of  pressure  on  volumetric  strain  and  the  depen¬ 
dence  of  deviatoric  stress  on  shear  strain,  the  pressure  will  also  generally 
depend  on  the  shear  strain,  and  the  deviatoric  stress  will  also  generally 
depend  on  the  volumetric  strain.  Furthermore,  these  latter  dependencies 
are  coupled.  Indeed,  from  (3.3)2t3,  (3.6)x ,  and  (3.7)x,  it  follows  that  the  co¬ 
efficients  p2(/5)  and  Pz{P)  in  the  expansion  (3.2)  for  the  pressure  are  related 
to  p  and  u  by 


Piifi)  =  P2 


d_ 

dp 


(3.9) 


U! 


(3.10) 


This  coupling  of  the  pressure  and  deviatoric  stress  may  be  made  more 
explicit  as  follows.  First,  note  that  from  (2.9)  and  either  (3.8)  or  (3.5),  we 
have 


IIT‘113  =  («,-)!  +  (t;y  +  (is?  =  [MPlYf  +  CMIMl3) .  (3.n) 

Then  on  solving  (3.11)  for  7*  and  substituting  into  (3.2)2,  we  obtain  the 
following  relation  between  the  pressure  p,  the  hydrostatic  pressure  p  =  p(p), 
the  shear  modulus  p,  and  the  equivalent  shear  stress  ||T*||: 

p  =  pW  +  4^(^-'‘)||t’|Ij+^(IIt'H3) 

=  ' + ^  (f  s  -  0 
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where  in  (3.12)x  the  coefficient  of  ||T*||2  is  regarded  as  a  function  of  p,  and 
in  (3.12)2  it  is  regarded  as  a  function  of  p. 

These  results  may  also  be  expressed  in  terms  of  the  deviatoric  parts 
E*,  G*,  and  H*  of  the  (Eulerian)  finite  strain  tensors 

E  =  V-I,  GsJ(B-I),  HsJ(I-B->);  (3.13) 

1 

in  particular,  H  is  known  as  the  Almansi-Hamel  strain  tensor.  Some  care 
must  be  taken  here  since,  unlike  L*,  the  tensors  E*  =  V*,  G*  =  |B*,  and 
H*  =  )*  are  not  independent  of  the  volumetric  strain.  Analogous 

to  the  definition  (2.14)  of  L**,  for  any  tensor  A  let  A**  =  [(A*)2]*  = 
(A*)2  —  |tr[(A*)2]I.  Then  it  can  be  shown  that 

L*  w  p/*Em  -  §p2/3E** 

»  p2/3G*  -  p4/3 G**  «  p~V3H*  +  p~4/3H**  ,  (3.14) 

L**  »  pV3 E~  »  p4l3G**  w  p~4/3 H** ,  (3.15) 

7*  «  p2/3||E*||2  »  /54/3||G*||2  «  p_4/3||H*||2  ,  (3.16) 

where  the  error  in  each  of  these  approximations  is  of  order  ||L*  ||3.  Substitu¬ 
tion  of  (3.14)-(3.16)  into  (3.1)x,  (3.2)2,  and  (3.8)2  yields  alternate  formulas 
for  e,  p,  and  T*  to  within  an  error  of  order  ||L*||3.  Also  note  that  if  the 
E**,  G**,  and  H**  terms  are  omitted  in  (3.14),  then  the  error  in  (3.14)  is 
of  order  ||L*||2,  and  (3.8)x  and  (3.2)x  yield 

T*  »  2p(p)L*  «  2pV3p(p)E* 

»  2p2'3ii(p)G*  «  2p~2!3 p(pya*  (3.17) 

and  p  p(p),  all  to  within  an  error  of  order  ||L*||2.  Thus  when  terms 
of  second  order  in  the  shear  strain  are  neglected,  the  pressure  depends 
only  on  the  density,  whereas  the  deviatoric  stress  depends  on  both  the 
shear  strain  and  the  density  through  the  density  dependence  of  the  shear 
modulus.8  As  we  see  in  the  next  section,  for  many  materials  the  shear 
modulus  changes  substantially  over  the  range  of  densities  encounterd  in 
high  velocity  impacts. 

8The  reader  who  prefers  to  work  with  the  strain  tensor  G,  for  example,  may  be 
tempted  to  use  the  term  “shear  modulus”  for  the  coefficient  p2l3p(p)  of  2G*  in  (3.17). 
However,  as  we  see  in  the  next  section,  this  would  be  in  conflict  with  the  terminology 
used  in  much  of  the  wave  propagation  literature. 
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4  Material  Properties 


Let  Ul  and  U$  denote  the  longitudinal  and  shear  wave  speeds  in  an  isotropic 
hyperelastic  material  in  a  state  of  dilatational  strain  (cf.  (2.13))  under 
hydrostatic  stress  T  =  —p(p) I.  These  wave  speeds  are  related  by  the  well- 
known  formula  (cf.  [1,  eq.  (75.4)]) 


where 


is  referred  to  as  the  bulk  wave  speed.  Then  from  (3.4)  we  have 


(4.1) 


(4.2) 


K  =  pU*.  (4.3) 

Similarly,  it  is  customary  to  define  the  shear  modulus  p  (regarded  as  a 
function  of  p ,  p,  or  p)  by 

p  =  p  Us2.  (4.4) 

Then  | p  =  p  —  p  dp/ dp,  and  by  use  of  the  well-known  formula  for  the 
longitudinal  wave  speed  (cf.  [1,  eq.  (74.7)]), 

M  =  W  («) 


it  can  be  verified  that  the  shear  modulus  as  defined  here  agrees  with  the 
definition  used  in  the  previous  section.  The  shear  and  bulk  moduli  as 
a  function  of  density  or  hydrostatic  pressure  are  usually  obtained  from 
measurements  of  ultrasonic  longitudinal  and  shear  wave  speeds,  together 
with  the  above  relations.  From  (4.4),  (3.9),  and  (3.4),  we  see  that  the 
coefficient  p2  of  7*  in  the  expansion  (3.2)  for  p  is  also  given  by 


p2  =  2pnUs 


dUs_ 

dp 


Then  from  (3.9)  and  (4.6)  it  follows  that 


dp 

p2>  0  <£>  p—  >  p  & 
dp 


dpL>P 

dp  k 


d/Us_ 

dp 


0. 


(4.6) 


(4.7) 


In  particular,  dUs/dp  >  0  implies  dp/ dp  >  0.  Also,  since  7*  >  0  for  any 
nonzero  shear  strain,  from  (4.7)  and  (3.2)2  we  see  that  if  dUs/dp  >  0,  then 
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shear  strain  at  constant  density  results  in  an  increase  in  pressure,  whereas 
shear  strain  at  zero  pressure  results  in  bulking  (since  p(p)  <  0  in  this  case), 
provided  that  terms  of  order  three  in  the  shear  strain  can  be  neglected.  For 
most  isotropic  elastic  materials,  the  shear  wave  speed  does  in  fact  increase 
with  pressure.9 

The  elastic  moduli  and  their  pressure  derivatives  evaluated  at  p  =  0 
(equivalently,  p  =  1)  are  denoted  with  a  zero  subscript.  Values  for  a  few 
materials  are  listed  in  Table  l.10 

Table  1:  Elastic  Moduli  and  Their  Pressure  Derivatives  (k0  and  p0  in  GPa, 
(dp/dp) 0/ p0  in  GPa-1,  other  quantities  dimensionless). 


Material 

K0 

Vo 

O 

se  |<a, 
Its 

o 

^3 

wmm 

It: 

o 

1 

*  1^ 

O  I  o 

Al 

76.0 

26.1 

4.4 

1.8 

0.069 

1.5 

Cu 

137 

47.7 

5.5 

1.4 

0.028 

1.0 

U 

113 

84.4 

5.9 

3.0 

0.036 

2.3 

W 

310 

160 

3.9 

2.3 

0.014 

1.8 

TiB2 

237 

246 

2.0 

9.0 

0.036 

7.9 

ZnO 

139 

44.2 

4.8 

-0.69 

-0.016 

-1.0 

fused  silica 

36.7 

31.3 

-6.3 

-3.2 

-0.10 

-4.0 

Table  1  shows  that  poly  crystalline  zinc  oxide  and  fused  silica  are  anoma¬ 
lous  in  that  their  shear  modulus  decreases  with  increasing  pressure.  For 
these  materials  the  inequalities  in  (4.7)  are  reversed,  in  which  case  (3.2) 
implies  shear  induced  compaction  (for  sufficiently  small  shear  strains)  at 
zero  pressure.  This  is  indeed  a  well-known  phenomenon  in  fused  silica. 
Observe  that 


dp 

Vo  +  ~TZ 
dp 


p  +  0(p2)  and 


V  ~/*o 
Vo 


1  dp 
Vo  dp 


p  +  0(p>).  (4.8) 


Thus  (dp/ dp)0f p0  is  a  measure  of  the  relative  change  in  the  shear  modulus 
due  to  hydrostatic  pressure. 


9This  statement  requires  some  qualification  in  the  adiabatic  case,  since  the  shear 
modulus  and  the  shear  wave  speed  should  increase  with  pressure  to  some  peak  values 
and  then  decay  to  zero  as  the  material  melts.  Also  note  that  (4.7)  remains  valid  if  all 
inequalities  are  replaced  by  equality,  in  which  case  p  =  p0p-,  thus  for  the  typical  case 
where  dUs/ dp  ^  0,  the  shear  modulus  (and  hence  T  )  cannot  be  proportional  p. 

10The  data  for  Al,  Cu,  U,  and  W  are  taken  from  [2],  TiB2  from  [3],  ZnO  from  [4], 
and  fused  silica  from  [5]. 
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Expansion  of  the  hydrostatic  pressure  p  in  powers  of  p  yields 
Pit)  =  K°(P  ~  V  +  Y  (%  ~  x)  (^  “  !)2  +  ^((^»  “  !)3)  •  (4-9) 

Now  consider  the  case  where  ||L*||  =  0(p  —  1);  i.e.,  where  the  shear  strain 
and  the  volumetric  strain  are  of  the  same  order.  Then 

7*  =  a  (p  —  l)2  +  G((p  -  l)3)  (4.10) 

for  some  constant  a  >  0  which  depends  on  the  particular  deformation  in 
question;  and  from  (3.2),  (3.9),  and  (4.9),  we  have 


Equation  (4.11)  gives  an  estimate  for  the  relative  increase  in  the  pressure 
due  to  shear  strain.  Equation  (4.12)  is  useful  for  comparing  the  second- 
order  contributions  to  the  total  pressure:  the  first  group  of  terms  inside 
the  square  brackets  is  the  contribution  due  volumetric  strain;  the  second 
group  is  the  contribution  due  to  shear  strain. 

These  results  apply  in  particular  to  uniaxial  strain  where  A2  =  A3  =  1, 
in  which  case  7*  =  |(ln  p)2  «  |(p— l)2,  and  hence  a  =  2/3  in  (4.10)-(4.12). 
We  will  apply  the  above  results  to  uniaxial  strain  of  the  polycrystalline  ti¬ 
tanium  diboride  tested  in  [3].  From  the  data  in  Table  1  we  find  that  the 
contributions  of  volumetric  strain  and  shear  strain  to  the  second-order  term 
in  the  expansion  (4.12)  for  the  total  pressure  are  0.5  and  5.27,  respectively; 
i.e.,  the  second-order  shear  strain  contribution  to  the  total  pressure  is  an 
order  of  magnitude  greater  than  the  second-order  volumetric  strain  con¬ 
tribution.  And  from  (4.11)  we  have  ( p  —  p)jp  fa  5.27  (p  —  1).  Of  course, 
these  results  are  valid  only  within  the  elastic  range  of  this  material.  At 
the  Hugoniot  elastic  limit  (HEL),  p  —  1  is  approximately  equal  to  0.03  (cf. 
[6]),  so  that  (p  —  p)/p  «  0.16.  In  other  words,  at  the  HEL  the  relative 
increase  in  pressure  due  to  shear  strain  is  about  16%.  If  not  taken  into 
account,  this  increase  would  lead  to  an  overestimation  of  the  shear  stress 
at  the  HEL  (obtained  by  subtracting  p  from  the  longitudinal  stress). 
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